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We are in motion against the cosmic backdrop. This motion is evidenced 
by the systematic temperature shift — or dipole anisotropy — observed in 
the Cosmic Microwave Background radiation (CMB) |IJ. Because of the 
Doppler effect, the temperature of the CMB is 0.1 per cent higher in our 
direction of motion through the Universe. If our standard cosmological 
i/-) ' understanding is correct, this dipole should also be present as an enhance- 

ment in the surface density of distant galaxies 0. The main obstacle in 
£T) ■ finding this signal is the very uneven distribution of nearby galaxies in 

the Local Supercluster, which drowns out the small cosmological imprint. 
Here we report the first detection of the expected dipole anisotropy in the 
galaxy distribution, in a survey of galaxies detected in radio waves. Radio 
galaxies are mostly located at cosmological distances, so the contamination 
from nearby clusters should be small. With local radio sources removed, 
we find a dipole anisotropy in the radio galaxy distribution in the same 
direction as the CMB, close to the expected amplitude. This result is 



confirmation of the standard cosmological interpretation of the CMB. 



The dipole anisotropy in the galaxy distribution due to our motion - the "velocity 
dipole" - has two causes. Firstly, emissions from galaxies ahead of us are Doppler 
boosted with respect to those behind, enhancing their fluxes; hence more are detected 
above a given flux threshold. Secondly, aberration of angles (as predicted by special 
relativity) causes the position vectors of galaxies to be displaced towards our direction 
of motion. We must distinguish the velocity dipole - an imprint in the distribution of 
distant, evenly distributed galaxies - from the anisotropic distribution of local galaxies 
due to the Local Supercluster. We describe the latter as the "clustering dipole"; it 
will not be a dipole distribution, but will have a dipole component. The clustering 
dipole is the cause of our motion through the Universe (the gravitational pull of the 
Local Supercluster), whilst the velocity dipole is the consequence of it. To measure 
the velocity dipole we must remove the contamination of the clustering dipole. 

The velocity dipole is a small effect, causing a surface density enhancement of ~ 
1 per cent. In order to detect it we must identify a sufficient quantity of distant 
objects and maintain uniformity in calibration to better than 1 per cent over the 
sky. Active galactic nuclei (AGN) are situated in galaxies at cosmological distances 
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and are therefore ideal probes of the velocity dipole. Unsuccessful attempts have 
been made to detect it in AGN X-ray surveys |J and radio surveys Q. In contrast, 
the clustering dipole has been detected and quantified in studies probing the local 
distribution of galaxies at both optical || and infra-red || wavelengths. In particular, 
Rowan-Robinson et al. used the Infra-Red Astronomical Satellite (IRAS) PSCz 
catalogue to compute the implied motion of the Local Group of galaxies. 

We measure the galaxy dipole using the NRAO VLA Sky Survey of radio sources 
(NVSS) 0. The NVSS covers a large fraction of the celestial sphere (82 per cent) 
and is deep (~ 50 sources deg -2 ), thus providing sufficient sources for good statistics. 
Moreover, the radio AGN in the survey are at a median redshift z ~ 1 ||, so that 
the contamination from local sources should be minimal. To prevent pollution from 
Galactic radio sources, we restrict our analysis to Galactic latitudes greater than 15°. 
The expected dipole amplitude is so small that systematic fluctuations above 1 per 
cent over wide angles will drown dipole signal. At flux-densities below 15 mJy, the 
NVSS does suffer from spurious systematic fluctuations in source surface density of 
this magnitude (Figure p. 

To eliminate the clustering dipole we remove from the sample radio sources within 
30 arc-seconds of known nearby galaxies as listed in the IRAS PSCz catalogue 
and the Third Reference Catalogue of Bright Galaxies (RCBG3) [It|. The IRAS 
galaxies eliminate local spiral galaxies obeying the far- infrared/radio correlation, and 
the RCBG3 sources remove radio sources powered by AGN in local elliptical galaxies. 

To measure the dipole from the remaining distribution of galaxies, we measure the 
spherical harmonics of the distribution. Any density field across the sky a(6, <fi) - 
where 9 (0 — > it) and <ft (0 — > 2tt) are ordinary spherical polar co-ordinates - can be 
expanded as a sum of spherical harmonics Yi tTn (9, <f>) with coefficients A^ m \ 

oo +Z 

<KM)=E E (i) 

1=0 m=-l 

The Ith spherical harmonics produce fluctuations on angular scales ~ 180°// and are 
described by 21 + 1 independent coefficients. A dipole distribution across the sky will 
produce signal in the three / = 1 spherical harmonic coefficients, but no signal in 
the others. To estimate the coefficients A^ m from N discrete objects with positions 
(9i,(pi), we calculate 

N 

A l , m = J2Y i : m (0 i , ( j>i) (2) 

i=l 

We measure the harmonic coefficients A^ m of the NVSS galaxy distribution for 1 < 
/ < 3 using equation || (a total of 15 independent coefficients). We measure the 
I > 1 coefficients because, as the sky is incomplete (containing blank regions below 
declination —40° and inside Galactic latitude 15°), higher harmonics are influenced 
by a dipole distribution. Further, any calibration gradients would (in general) bias all 
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harmonics. Having measured the coefficients, we then determine how well they are fit 
by a dipole model (with the same blank regions). 



A dipole model has three parameters: an amplitude 5 and direction described by 
two angles (#,</>). We define 5 as the total fractional fluctuation in surface density 
from maximum to minimum. Table [I] displays the best-fitting dipole parameters and 
errors for a sequence of flux-density thresholds and describes the method used. The 
dipole model is a good fit (Xred ~ 1) down to 15 mJy Below this level the NVSS 
data- reduction biases of Figure [I] become significant, producing a sudden increase in 
Xred an d a significant shift in the best-fitting 9 (but not 0), anticipated because the 
biases are declination- dependent. Figure |2] illustrates the best-fitting direction (9, (j>) 
of the measured dipole at intensity levels above 20 mJy. The direction does not change 
with flux-density threshold and is in good agreement with the direction of the CMB 
velocity dipole 9 = 97.2 ± 0.1°, <p = 168.0 ± 0.1° |T[J. The local clustering dipole 
typically contributes A5 = 0.5 x 10~ 2 to the total amplitude. 

Ellis & Baldwin 0] derived the expected amplitude of the velocity dipole assuming 
radio sources to have identical power-law spectra S(u) oc v~ a and an integral source- 
count N(> S) oc S~ x : 

5 pred = 2(v/c)[2 + x(l + a)] (3) 

where v is the peculiar velocity of our reference frame with respect to the frame in 
which the radio source population is (assumed) isotropic, and c is the speed of light. 



CMB dipole measurements imply v = 370 ± 2 km s Jl ![]. Taking x ~ 1 and a mean 
spectral index a ~ 0.75, we find S pre d = 0.9 x 10~ 2 . Note that v is the peculiar velocity 
of the Sun in the CMB frame. (The speed of the Earth around the Sun is ~ 30 km 
s _1 and can be neglected). Our observed values of 5 are on average 1.5 a away from 
this prediction. 

We can be confident that our measured surface- density dipole is due to velocity rather 
than residual local clustering. Almost all the contribution to the local clustering dipole 
comes from redshifts z < 0.03 - analysis of the IRAS PSCz dipole || shows that the 
dipole amplitude has almost converged by 100 h^ 1 Mpc. This result is verified in 
Figure |3|, where we plot the strength of the dipole contribution in redshift shells from 
the PSCz and RCBG3 catalogues. Furthermore, these catalogues contain almost 
all the radio galaxies to z — 0.03. This is demonstrated by Figure |] which gives 
magnitude-redshift plots for PSCz and RCBG3 galaxies matched with NVSS sources. 
The distributions make it clear that we are missing a negligible fraction of radio 
galaxies at z < 0.03. 

It is hard to conceive of any other effect that can produce our measured dipole. A 
signal originating from instrumental or data reduction effects would tend to skew the 
dipole direction towards the north or south poles, as such effects only depend on 
declination. However, our measured dipole points to declination pa 0° in agreement 
with the CMB dipole. Moreover instrumental effects would (in general) produce 
signal in higher (non-dipole) harmonics, but the best-fit dipole model is a good fit 
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with xled ~ 1- Furthermore, the dipole amplitude is extremely unlikely to be a 
chance fluctuation of an isotropic universe. For example, for a flux-density threshold 
of 20 mJy, the isotropic universe (5 = 0) model fits the observed harmonics with 
Xred = 2.28, rejecting this model at a confidence level of > 99.5 per cent. It is difficult 
to see what else could produce the dipole, particularly one pointing (within the la 
direction contours) in the CMB dipole direction. For the 20 mJy measurement, the 
la confidence region for the direction encompasses just 5 per cent of the sky. 

The alignment of our measured velocity dipole with the CMB dipole supports the 
predictions of fundamental assumptions about the cosmological origin of the CMB, 
the lynchpin of much present-day observational cosmology. This result indicates that 
objects at z ~ 1 (such as distant radio galaxies) can define the cosmic frame and have 
an isotropic distribution in that frame. 
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Table 1: 



Flux 


N 


Best S 


Best 


Best 9 


"X-red 


/ mJy 




xlO" 2 


1° 


1° 




> 40 


125,603 


1.4 ± 0.9 


149 ± 49 


135 ± 38 


1.02 


> 35 


143,524 


1.7 ± 0.7 


161 ± 44 


117 ± 39 


0.74 


> 30 


166,694 


2.2 ± 0.7 


156 ± 32 


88 ± 33 


1.01 


> 25 


197,998 


2.2 ± 0.6 


158 ± 30 


94 ± 34 


1.01 


> 20 


242,710 


2.1 ± 0.6 


153 ± 27 


93 ± 29 


1.32 


> 15 


311,037 


1.5 ± 0.5 


148 ± 29 


59 ± 31 


1.81 


> 10 


431,990 


1.0 ± 0.4 


132 ± 29 


25 ± 19 


4.96 



Table |T|: The dipole amplitude and direction that best fit the galaxy dis- 
tribution. The Table lists the best-fitting dipole model parameters (5, 8, eft) together 
with la errors for NVSS sources above various flux thresholds, after local sources 
have been removed. 5 is expressed as a percentage. N is the number of sources left 
in the unmasked survey region and the reduced \ 2 values of the fit to the observed 
survey harmonics are also shown. The expected (CMB) velocity dipole parameters 
are S = 0.9 x 10~ 2 , <p — 168°, 9 = 97°. To find the best-fitting parameters, we vary 
the values of (S, 6, (ft) over a grid. For each set of parameters, we randomly generate a 
dipole galaxy distribution and mask it in the same way as the real data (below decli- 
nation —40° and inside Galactic latitude 15°). We then measure the Ai^s of this test 
distribution, averaging over 100 realizations to increase accuracy. We hence obtain a 
mean and standard deviation for each A^ m in the set corresponding to the parameters 
(8, 9, (ft). We can then compute the x 2 statistic with the observed survey harmonics. 
The best-fitting dipole parameters are found by minimizing x 2 ; this minimum value 
(goodness-of-fit) indicates how well the dipole model fits the data. The errors on 
the best-fitting parameters are properly described by x 2 contours in the parameter 
space. For example, with all three parameters varying, the la and 2a error regions 
are enclosed by \ 2 increases of 3.53 and 8.02. 
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Figure [I]: The NVSS suffers from systematic fluctuations in source surface 
density across the sky at low flux-densities. Here we plot the surface density as 
a function of declination, in declination bins of width 10°, for flux-density thresholds 
3.5 mJy (filled circles) and 15 mJy (open circles). The declination range of each 
array configuration (D or DnC) is shown. The 2 per cent fluctuations present at 
the completeness limit 3.5 mJy largely disappear when the threshold is raised to 15 
mJy. These biases result from the sparse «i>-plane sampling of the NVSS (W. Cotton, 
private communication, 2001). The fluctuations are a function of declination only, as 
the NVSS beam and analysis procedure have no way of knowing the right ascension 
of objects. The NVSS was carried out with the Very Large Array at 1.4 GHz in D 
and DnC configurations during 1993 to 1996. The survey covers all the sky north 
of declination —40°; the catalogue contains ~ 1.8 x 10 6 sources and is claimed to be 
~ 99 per cent complete at integrated flux density 3.5 mJy. The Full-Width Half- 
Maximum of the synthesized beam is ~ 45 arcsec. Prior to analysis, we masked 22 
regions around bright and extended radio galaxies; these can appear in the NVSS 
catalogue as clusters of many point sources by virtue of the source-finding algorithm. 
These regions constitute a tiny fraction of the surveyed area and were filled with 
random distributions of sources at the mean surface density of the rest of the survey. 
The large NVSS beam means that only a small fraction of galaxies are resolved into 
double sources - this is important because a large number of doubles would increase 
the source surface density and hence the observed dipole amplitude. 

Figure ^: Our measured dipole direction agrees with that of the CMB 
dipole. Here we plot la and 2a contours for the best-fitting dipole direction for flux- 
density thresholds of 40 mJy (dotted), 30 mJy (dashed) and 20 mJy (solid). These 
contours are defined by A% 2 = 2.30 and 6.17 (with 6 and <fi varying and 5 fixed at its 
best-fitting value). The filled circle illustrates the direction of the CMB dipole, which 
is consistent with our measurements. 

Figure |3|: The clustering dipole originates from redshifts z < 0.03. Here we 
measure the strength of the clustering dipole in the PSCz and RCBG3 catalogues as 
a function of redshift. At each redshift z, we average over a redshift range Az = 0.01 
and thus the plot starts at z = 0.005. For each redshift shell, the y-axis plots v ' D 2 — N 
where D = for the galaxies in that shell [|J. This expression is proportional 

to the imbalance in the total number of sources across the sky, the quantity to which 
the total dipole amplitude is sensitive. The peak in RCBG3 at z — 0.025 is due to 
the Coma cluster. Hence it seems that once we get past Coma, the local clustering 
dipole is negligible. 

Figure ^|: The PSCz and RCBG3 catalogues include almost all nearby 
NVSS galaxies. Here we present magnitude-redshift scatter plots for PSCz galaxies 
(left-hand panel) and RCBG3 galaxies (right-hand panel) matched with NVSS sources 
brighter than 15 mJy. The horizontal lines are the stated completeness limits for the 
two catalogues. Clearly the bulk of the radio galaxies lie safely above the completeness 
limits up to at least z = 0.03. Although the number of RCBG3 matches starts to die 
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away at z 0.025, this is likely due to the N(z) of the radio population rather than 
any serious incompleteness in RCBG3 - doing the same plot for all RCBG3 galaxies 
shows no such fall-off until the magnitude threshold is reached. 
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Figure 1: 
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Figure 2: 
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Figure 4: 
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